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We re-examine the thermal evolution of the postshock layer in primordial gas clouds. 
Comparing the time scales, we find that the evolutionary paths of postshock regions in 
primordial gas clouds can be basically understood in terms of the diagram drawn in the ion- 
ization degree vs temperature plane. The results obtained from the diagram are independent 
of the density in the case that we do not include photodissociation and photoionization. We 
also argue that the diagram is not only relevant to the case of the steady postshock flow, 
but also to the isochorically cooling gas. 

§1. Introduction 

The thermal evolution of primordial gas clouds has been investigated by many 
authors. Eif Almost all of those studies have been concerned with 

the formation of various kinds of galaxies, or primordial stars. In those papers, 
galaxies are assumed to grow out of small density perturbations present in the early 
universe. Because of the coldness of the growing density perturbations, the formed 
clouds, which are the progenitors of galaxies, experience strong shock heating at 
the bouncing epoch. Shock heating is also expected in the hierarchical clustering 
scenario of structure formation. In this case, shocks are expected in the collision 
between two clouds which are trapped in the gravitational potential well associated 
with the larger structure. 

In any case, shock heating is expected in the era of galaxy formation. The spa- 
tial structure and the thermal evolution of the postshock layer in primordial gas 
clouds were investigated by many authors. I' ©®S>H>0)>&>§> In those stud- 
ies the common and the most important pointis the over production of hydrogen 
molecules. For example, in Shapiro and Ranged (hereafter SK), the thermal evo- 
lution of steady postshock flow is investigated. They found that the postshock flow 
cools down so fast that the recombination process cannot catch up with the cooling. 
As a result, the ionization degree remains high (y e ~ 10 -3 ), even when the temper- 
ature has dropped below 10 4 K. Feeded these "relic" electrons, hydrogen molecules 
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form through the processes 

H + e^H"+7, (1-1) 
H + H"^H 2 + e. (1-2) 

SK stressed the importance of these non-equilibrium cooling processes and the dif- 
ficulties of describing these processes with a finite-difference scheme. 

SK calculated the cases corresponding to the shock velocity in the range (50km/s— 
300km/s). In all cases, the resulting fractions of hydrogen molecules are ~ 10~ 3 , ir- 
respective of the initial shock velocities. Ionization degrees are also independent of 
the initial shock velocities for the same temperature if T <C 10 4 K. Surprisingly, the 
convergence is quite good, though chemical equilibrium is not achieved throughout 
almost the entire evolutionary path. Therefore there may be a clear explanation for 
the independence of the initial shock velocity. 

In this paper we re-examine the thermal evolution of the postshock layer in pri- 
mordial gas clouds. To elucidate the convergence of the hydrogen molecule fraction 
in the course of the evolution of the postshock region, we simply compare several 
important time scales (§|2|). We find that we can predict the thermal evolution of 
the postshock layer for any initial conditions in a diagram of the ionization degree 
versus temperature. In §|| we also test the results obtained in §|| by comparing the 
results with those of numerical calculations for the steady postshock flow. In §|| we 
show that the thermal evolution of isochorically cooling gas is very similar to that 
of steady postshock flow. The remainder of the section is devoted to the further 
discussion and summary. 

§2. TIME SCALES 

In this section we compare the time scales which are essential to the thermal 
evolution of the postshock layer in primordial gas clouds. First of all, we enumerate 
the definition of the time scales 

pkT 
' iinipA ' 

1 



^cool — Cp jV A , (2T) 



k ion n H n e k ion n N (l-y e ) 



(2-2) 



_ n p 1 
k TCC n e n p k vcc n^y e 

tdis = £§3 = * , (24) 

Eij Hf "<".; Eij kjf "Ml HI. i 

Here t coo \, t lon , t Tec , tdis and tf or denote the time scales of cooling, ionization, recom- 
bination, H2 dissociation, and H2 formation, respectively, til tin, n-i, yi (= nj/njv), 

*^ The time scale of thermalization in the postshock layer is shorter than the other time scales. 
For example, at 2 x 10 4 K the time scales of thermalization for plasma and neutral gas are about 
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k x , p, p and m p are the number density of nucleons, the number density of the zth 
species, the fraction of ith species, the chemical reaction rate coefficient of the "X" 
process, the mass density, the mean molecular weight, and proton mass, respectively. 
A denotes the net cooling rate which includes radiative and chemical cooling/heating, 
assuming that the system is optically thin for cooling photons. This assumption is 
valid for subgalactic clouds with mass 10 6 — 1O 8 M0 and whose densities are smaller 
than ~ 10 8 cm -3 . If a cloud violates this condition, the system could be optically 
thick for the H2 cooling line photons. Another assumption is that the primordial 
gas does not contain helium but only hydrogen. This assumption affects the mean 
molecular weight and the cooling rate at high temperature (T ~ 10 5 K). However, 
these are minor effects for the thermal evolution of primordial gas clouds at lower 
temperature (T Ss 10 4 K), which we are especially interested in. Therefore, for sim- 
plicity we do not include helium. C PjV denotes the heat capacity of the gas cloud. If 
the gas evolves isobarically, C p is chosen, and if the cloud is isochoric, C v is chosen. 
These quantities are expressed in terms of 7 as 

C P = ^ (2-6) 

C v = (2-7) 

where 7 = 5/3 for an ideal gas without molecules. When we are interested in a 
steady flow, the postshock layer is almost isobaric (see Appendix |B|). 

The ionization time scale (Eq. (2-2)) is defined as the variation time scale of y e . 
Another possibility is definition in terms of 1 — y e . However, it is not appropriate for 
the problem we are now interested in, because the magnitude of ionization degree 
itself is important for the formation of hydrogen molecules. In any case, comparing 
these five time scales we can find the fastest process without solving detailed time 
dependent differential equations. 

We should remark that all of these time scales are proportional to p^ 1 if the 
cooling rate and other reaction rates are proportional to p 2 . The dominant compo- 
nent of the cooling changes according to the temperature. Below 10 4 K, the H2 line 
emission dominates the total cooling of the cloud. In this process, the cooling rate is 
proportional to the number density. oLthe hydrogen molecules at the excited levels 
whose fraction is proportional to p.OP til Hence, the cooling time scale is proportional 
to p- 1 for T & 10 4 K as in Eq. (2-1). The cooling process in 10 4 K ^ T ^ 10 5 K 
is dominated by the bound-bound transition of the hydrogen atoms, and the num- 
ber of the excited hydrogen atoms is determined by the collisions with other atoms. 
Therefore, the cooling rate is proportional to the square of the total density. For 
higher temperature (T <; 10 5 K), the free- free emission by the collision between ions 
and electrons dominates the energy radiated away from the cloud. In any case, t^he, 
cooling rate is proportional to p 2 , and the cooling time is proportional to p~ 



10 2 and 10 5 times shorter than the ionization time scale, respectively. 

*' This treatment is limited to low density (n ^ 10 4 cm -3 ). The cooling rate is not proportional 
to p 2 for higher density (higher than the critical density). 

**' Compton cooling is effective at high redshift and high temperature. In this case, the cooling 
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The chemical reaction rates are proportional to p 2 if the reactions are dominated 
by the collisional processes. Photoionization and H2 photodissociation processes by 
UV photons emitted from the postshock hot region could affect the ionization rates 
and the H2 dissociation rates in the postshock regions. If these radiative reactions 
dominate the other reactions, total reaction rates will not be proportional to p 2 , be- 
cause photoionization and H2 photodissociation processes are not collisional. Kang 
and Shapiro E£P investigated this problem. The figures in their paper indicate that 
the effects caused by the internally produced photons are not too large to dominate 
that of the collisional processes. Hence, as a Oth-order approximation, we neglect 
the photoionization and the H2 photodissociation processes in this paper. 

As a result, the ratio of any two time scales is independent of the total density, 
and they are determined by the temperature and the chemical compositions (y e 
and j/h 2 )- As will be discussed later, yu 2 is almost always determined by y e and 
T. Therefore, the ratio of any two time scales just depends on T and y e . In the 
following sections, we compare the time scales individually. 

2.1. Ionization and cooling 

First, we investigate the case ti on < ^reo which corresponds to higher tempera- 
ture (T ^ 10 4 K), and we compare the time scale of ionization with that of cooling. 
As mentioned previously, the ratio ti on /t coo i does not depend on p, but on the ion- 
ization degree y e and temperature T, because both ionization and cooling, which 
are dominated by atomic hydrogen line cooling, occur through two-body collisions. 
Hence the ratio iionAcool is determined at each point on the y e -T plane indepen- 
dently of the density. The line along which t\ oa /t coo \ is equal to unity is indicated by 
the solid line in Fig. 1. If the initial conditions are given below the line, the system 
is ionized before it cools down, because t\ on /t coo \ < 1. Otherwise, the system cools 
before it is ionized. The expected evolutionary path is also expressed by the arrows 
in Fig. 1. 

2.2. Recombination and cooling 

Next we investigate the case t rec < t; on , which corresponds to lower tempera- 
ture, and examine the balance between the recombination process and the cooling 
process. The recombination process is important for T <, 10 4 K . When T <, 10 4 K, 
the line cooling of the hydrogen molecules dominates the atomic hydrogen line cool- 
ing. Therefore the ratio t rec /t coo \ depends not only on the ionization degree and the 
temperature but also on the fraction of hydrogen molecules. In order to draw the 
line t Tec /t coo \ = 1 in the y e -T plane, we need information on the fraction of hydrogen 
molecules. As will be discussed in the next subsection, the fraction of H2 is deter- 
mined by the chemical equilibrium above ~ 10 4 K. Then, it is approximately frozen 
below ~ 4 x 10 3 K through the evolutionary course with cooling. The line along 
which trecAcool = 1 is satisfied is drawn in Fig. 1. In the region above the line, the 
system recombines before it cools down, and below the line, the system cools down 
before the recombination proceeds. The expected evolutionary path is also indicated 
in Fig. 1 by the arrows. 



rate is not proportional to p 2 , but this does not change the results significantly. 
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Fig. 1. The y e -T plane is divided into several regions in which the ratio of the time scales are differ. 
The arrows in the plane denote the directions of the evolution of the postshock layer. The dot- 
short-dashed curve where the equality t lac = tion holds. Region (a) in which the condition 
iion < tcooi holds is bounded by the solid line. In this region, the system evolves upward in 
the y e -T diagram. In region (b), the conditions t coo \ < ti on , t lac always holds. Therefore the 
temperature drops down before the system is further ionized or recombines. In region (c) we 
have tree < tcooi, and the system evolves downward. Region (c) is bounded by the dotted line, 
short-dashed line and the long-dashed line. Those three lines denote the boundary on which 
the condition t coo i = t rcc holds. The long-dashed line was calculated by assuming that the 
cooling rate is dominated by the hydrogen line cooling. The dotted line was calculated assuming 
chemical equilibrium of the hydrogen molecules, and the short-dashed line was calculated with 
j/h 2 — 10 -3 , which corresponds to the frozen value of the hydrogen molecules for a shock velocity 
higher than ~ 30 km/s, see also Fig. 2. 

2.3. Fraction of H2 

Here we investigate whether or not H2 is in chemical equilibrium. If H2 is in 
chemical equilibrium (with given T and y e ), the cooling time is always determined 
by just T and y e . As a result* the thermal evolution of the system is completely 
determined in the y e -T plane.tJ The fraction of hydrogen molecules is roughly de- 
termined by estimating t^ s and tf or , which denote the time scale of the dissociation 
and the formation of hydrogen molecules. If at least one of the two time scales is 
shorter than the other time scales, t coo \ and t Tec , the fraction of hydrogen molecules 
changes toward the chemical equilibrium value. Hence, the ratios tdisAcool rec an d 
iforAcool,rec should be examined. 

First of all, we comment on the important properties of idis an d tf 0T . The time 
scale of dissociation is independent of yn 2 > t ne fraction of hydrogen molecules. On 

*' The other components such as H^" or H~ are always in chemical equilibrium, because they 
are very fragile. 
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the other hand, tf or is proportional to yu 2 . Then both the time scales given in terms 
of their equilibrium values are 

*dis = * d is' ( 2-8 ) 

ttor = &, (2-9) 

where the superscript eq denotes the value at which the hydrogen molecules are in 
chemical equilibrium for a given electron abundance and temperature. Here, the 
electron abundance is not an equilibrium value in general. Therefore the relation 
C = Cs holds. 

Suppose that we know the important time scale (t sys ) which characterize the 
evolution of the system. It could be trec^cooi, and so on. If t sys is independent of 
yu_ 2 i w e have general arguments as follows: 1) In the case yn 2 > Un eq i dissociation 
proceeds faster than the formation process. Hence, tdis controls the change of yn 2 . 
This is the reason we do not have to be concerned about tf or . Since t^is is independent 
of y H2 , y R2 converges to y^ if t dis < t sys is satisfied. 2) If y H2 < y^, y H2 is controlled 
by the time scale tf or , because the system will try to recover the fraction of H2 which 
is smaller than the equilibrium value. The maximal value of tf or is the equilibrium 
value, since tf or is proportional to yu 2 (Eq.(2-9)). Therefore, if t sys > tfo r (= tdis), ^sys 
is always greater than tf or . As a result, t sys > tdis is a sufficient condition for the 
system to maintain the chemical equilibrium of hydrogen molecules. In the following, 
we look into the concrete cases that the interesting time scales are t rec and t coo i- 

2.3.1. tdis, for VS t rec 

First, let us investigate the chemical equilibrium of hydrogen molecules compared 
to the recombination process. In other words, we investigate the magnitudes of 
^disArec and tfor/^rec- As discussed above, when we compare tdis and ty or with other 
time scales which do not depend on yn 2 , we do not have to be concerned about tf or . 
In this case, t rec is independent of yn 2 , and hence we only have to examine tdisArec- 
The time scale of recombination, t rcc (Eq.(2-3)), is proportional to y" 1 , and the time 
scale of dissociation is also proportional to y" 1 , because it is governed by collisions 
between H + and H2 . This means that the ratio tdis/^rcc 

is independent of y e . It is also 
independent of yn 2 , and hence the ratio just depends on the temperature. Therefore 
we have the critical temperature below which tdisArec is greater than unity (see Fig. 
2). This critical temperature is ~ 2600 K. This argument is different if the ionization 
degree is sufficiently low for the H — H2 collision to dominate the H + — H2 collision. 
However, within the regions we are interested (y e ^ 10 -5 — 10~ 6 ), the ionization 
degree is higher than this. The fact that the line t rec = tdis is almost vertical in Fig. 
2 confirms the validity of this treatment. 

2.3.2. tdis, for VS t CO ol 

Next, we study the chemical equilibrium of H2 with respect to the cooling of 
the gas; i.e. we consider the ratio tdi s /i C ooi and tforAcooi- However, we cannot use 
the previous argument in this case, because the cooling time scale depends on yn 2 
for T ^ 10 4 K. Accordingly, we must investigate the individual case again. 1) 
If yn 2 < y^ , we obtain the relation tf or < t^is- I n this case, we do not have to 
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Fig. 2. (a) The y e -T plane is divided into two regions, where chemical equilibrium is achieved for 
the fraction of hydrogen molecules and where not. The solid boundary of the shaded region is 
determined by the condition i coo i = £dis- The long-dashed line denotes the boundary at which 
tree = tdis is satisfied. In the shaded region, the chemical equilibrium of hydrogen molecules 
breaks down, and the fraction j/h 2 = nn 2 /nn freezes as the temperature and/or ionization 
degree drops. The dotted lines labeled as 100 km/s and 20 km/s denote the evolutionary path 
of y e . (b) The equilibrium and non-equilibrium fraction of hydrogen molecules are drawn in 
yn 2 -T plane. The solid lines labeled 100 km/s(j/H 2 ) and 20 km/s(j/H 2 ) denote the dynamically 
calculated evolutionary path of the fraction of hydrogen molecules. The short-dashed lines which 
are labeled as 100 km/s(y„ 2 ) and 20 km/s(j/g 9 2 ) denote the fraction of H2 calculated assuming 
the chemical equilibrium of hydrogen molecules with dynamically calculated electron density. 



consider t^is, because H2 has a smaller fraction than its equilibrium value. The ratio 
ifor/icool is almost proportional to y^ 2 , because t{ or is proportional to yu 2 and t coo \ 
is proportional to y^\ Therefore the ratio reaches a maximal value at yn 2 = y^ 2 . 
This means U ox /t coo \ is always smaller than t^/t e ^ oo[ . Then, if t^/t^, is smaller 
than unity, H2 is in equilibrium while the temperature drops. 2) If yu 2 > y^ 2 , we 
obtain tf or > tdis- In this case we just check the ratio idisAcooh because H2 has 
larger fraction than its equilibrium value. The ratio tdisAcool is proportional to ya 2 . 
Therefore, the ratio does not reach a maximum at 2/h 2 = 2/h 2 > ra ther a minimum. 
Then we cannot conclude that t^/t^^ < 1 is a sufficient condition for the chemical 
equilibrium according to the cooling process. However, in the region yn 2 > y^, the 
fraction of H2 converges to as the system cools, because the inclination of the 
short-dashed lines in Fig. 2(b) (denoting the evolution of y^ 2 ) is negative. If the 
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inclination is positive, yu 2 will depart from the equilibrium value as the temperature 
drops. In other words, if the hydrogen molecules are not in chemical equilibrium in 
the region j/h 2 > yn 2 i the cooling process proceeds until the system cools down to 
the state yn 2 = y^ 2 - 

In any case, if the condition idisAcooi < 1 is satisfied, the relation yn 2 — Vn 2 
always holds. In Fig. 2(a) the line on which the condition tdisAcool = 1 is satisfied is 
presented in the y e -T plane. The shaded region denotes the region in which chemical 
equilibrium for hydrogen molecules with respect to the cooling or recombination 
process does not hold. When the system comes into this region, the fraction of 
hydrogen molecules "freeze" . The fraction of hydrogen molecule becomes ~ 10~ 3 if 
the initial temperature is higher than ~ 2 x 10 4 K. This temperature corresponds 
to the initial shock velocity, ~ 30 km/s, when we discuss the steady postshock flow. 
For smaller shock velocities, the initial condition for the chemical compositions is 
not reset, and the fraction of hydrogen molecule becomes a smaller value. 

2.4. Shock diagram 

In the end of this section, we present a brief summary of §^. The main results 
of this section are summarized in Fig. 1. The evolution of a shock heated system 
is basically explained in the y e -T plane, which we call the "shock diagram" (Fig. 1). 
The evolutionary path of the postshock layer is obtained by tracing the directions 
of the arrows in the shock diagram. Postshock gas, which is heated up to high 
temperature, appears in region (a). In region (a), the shortest time scale is ii on and 
the system evolves upward and enters region (b). In region (b), i coo i is the shortest, 
and therefore the system evolves leftward this time and comes into region (c). In 
region (c), t rec is the shortest, and then the system evolves downward, and reenters 
region (b), where the system evolves leftward again and crosses the boundary of 
regions (b) and (c). Finally, the system evolves nearly along the short-dashed curve, 
which is the boundary of regions (b) and (c). For postshock gas which is not heated 
up to above 2 x 10 K, the thermal evolution is somewhat different from the previous 
case. In this case, when the system enters region (b) from region (a), the ionization 
degree is so low that the system does not enter region (c) at T ~ 10 4 K. Then the 
ionization degree and the fraction of hydrogen molecules do not "forget" the initial 
conditions of the chemical composition. 

2.5. Convergence of the ionization degree in the steady postshock flow 

The convergences of the ionization degree and the fraction of hydrogen molecules 
in the steady postshock flow are basically explained by Figs.l~3. At T ~ 10 4 K, the 
line t rec = icooi becomes nearly vertical, since the dominant cooling process is the 
hydrogen atomic line cooling which decreases by many orders below T~2x 10 4 K. 
If (1) the initial shock velocity is larger than ~ 30 km/s or (2) the initial ionization 
degree is larger than ~ 7 x 10 -3 and the initial temperature is larger than T ~ 10 4 K, 
then the evolutionary path should hit this nearly vertical line (the long-dashed line 
in Fig. 1). Because the line t Tec = t coo \ is nearly vertical at T ~ 10 4 K and y e ^ 10~ 2 , 
y e drops to 10~ 3 — 10~ 2 without changing the temperature and "forgets" the initial 
conditions, though y e is not in chemical equilibrium. The value converged to {y e ~ 
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Fig. 3. The evolution of the postshock layer in the y e -T plane is presented. The solid lines are 
the boundaries beyond which the direction of the evolution is changed, and each line has the 
same meaning as in Fig. 1. The short-dashed lines are the numerical results from the SK-type 
calculation. Each suffix denotes the initial shock velocity. 

7 x 10~ 3 ) is essentially determined by the equation t Tec (y e ,T) = i^ ol (y e ,T) at T ~ 
8000 K. Here, the temperature (8000 K) is obtained by the equality t ICC /t^^ = 1, 
where t Jf°[ is the cooling time estimated by the hydrogen line cooling. The equation 
^recAcooi = ^ depends both on y e and on T in general, but it depends just on T for 
Ve < 1. 

The convergence of H2 is also explained by Fig. 2. According to Fig. 2, H2 is still 
in chemical equilibrium just below T = 10 4 K, at which temperature the convergence 
of y e takes place. As a result, systems which satisfy the previous condition (1) or 
(2) experience the same state (same y e , same yn 2 ) just below T = 10 4 K. Therefore 
yu 2 also converges for different initial shock velocities. 

§3. COMPARISON WITH NUMERICAL CALCULATIONS 

In this section we compare the results obtained in the previous section with 
numerical calculations. The numerical calculations which we performed are the 
same as that in SK, except that our calculations do not contain helium and the 
cooling rate due to hydrogen molecules is based on Hollenbach and McKeeQM (see 
Appendices [A| and ||). 

Figure3 displays the evolution of y e as a function of the temperature behind the 



*' SK used the cooling rate in Lepp and Shull,^ but their formulae overestimate the cooling 
rate for lower temperature, T ^ 10 K. 
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shock front in steady flow. It is obvious that the expected evolutionary path in Fig. 
1 agrees very well with the numerically calculated path. We also find that there 
exists a critical value of the initial shock velocity above which the evolution of y e 
below 10 4 K converges. The critical value of the shock velocity is ~ 30 km/s. This 
condition for the initial shock velocity is equivalent to condition (1) in the previous 
subsection. 




Fig. 4. Evolution of the time scales in the steady flow calculation are presented. The initial shock 
velocity is f00 km/s. The solid line, dotted line, short-dashed line and long-dashed line represent 
tfon idis, ircc and t coo i, respectively. Every time scale is multiplied by density. 



The numerically calculated evolutionary path of t/h 2 is presented in Fig. 2. The 
chemical equilibrium for the hydrogen molecule holds outside of the shaded region 
in the y e -T plane. Once the system comes into the shaded region, the fraction of 
H2 departs from the chemical equilibrium and eventually freezes to a certain value, 
because the cooling proceeds faster than the dissociation and formation of hydrogen 
molecules. 

The evolution of the time scales in numerical calculation (v s = 100 km/s) is 
presented in Fig. 4. The dissociation time scale is identical to the formation time 
scale for T ^ 4000K, and they are smaller than the cooling time. For T 5s 4000K, 
tdis becomes larger than i C ooi> and yu 2 starts to depart from y# . At a slightly 
lower temperature, if or also becomes larger than t coo \. As a result, the value of yu 2 
approximately freeze for T 5s 4000K. The final value of yu 2 is 10~ 3 for v s ^ 30 km/s, 
and smaller for slower shock velocities. The frozen values degenerate for faster shock 
velocities, because the ionization degrees also degenerate for v s <; 30 km/s. 
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§4. Conclusion and discussion 

In this paper, we investigated the thermal evolution of the postshock layer in 
primordial gas clouds. Because of the delay of the recombination process compared 
to the line cooling due to the hydrogen molecules, significant amounts of electrons 
(y e ~ 10~ 3 ) remain below 8000 K. Feeded the electrons, H2 also maintains high 
abundance (yn 2 ~ 10 -3 ) for T 5s several x 10 3 K, then the cloud temperature drops 
to T ~ 100 K quickly due to the H2 line cooling. We have shown that both of the 
residual values of y e and yn 2 are universal for v s ^ 30 km/s by simple compari- 
son of time scales. The universality of the values to which y e and yn 2 converge is 
explained as follows: The evolution of the system is basically determined by T, y e 
and yn 2 . However, at high temperature, ya 2 attains its equilibrium value, and at 
low temperature, it will freeze to a value, because cooling precedes the dissociation 
and the formation of hydrogen molecules. Therefore the evolution of the system is 
predicted just by T and y e . The convergence of y e and yn 2 occurs when the coolant 
changes from hydrogen atoms to hydrogen molecules. In fact, assuming an isobaric 
cooling time, the predicted evolutionary paths agree very well with the numerically 
integrated paths for steady flow. 

These results are not restricted to the case of steady postshock flow, and what 
is more, the shock diagram describes not only the evolution of the postshock flow 
but also the primordial gas cloud which is once heated up to 10 4 K. The argument 
in §|^ simply depends on the comparison of time scales. When we compare the shock 
diagram (Fig. 1) with the numerically integrated results, the steadiness of the flow 
is used just as the isobaricity of the postshock flow in Eq. (2-1), by choosing C p . 
However, even if we choose C v in Eq. (2-1), i coo i is not much different from the 
cooling time for the isobaric flow. We have 

j.isobaric 

" C p /C v , (4-1) 



j-isochoric 
''cool 



= 7- (4-2) 

Consequently, the y e -T diagram obtained for the isochoric case, is very similar to 
Fig. 1. The diagram for the isochoric case is presented in Fig. 5. It is obvious that 
the difference between Fig. 1 and Fig. 5 is small. In the isochoric case, the final 
fraction of hydrogen molecules is ~ 10~ 3 again, as in the isobaric case. Hence, the 
intermediate case should have a similar diagram. Furthermore, the shock diagrams 
are applied not only to the postshock flow but also to primordial gas clouds which 
are initially sufficiently heated (T <; 10 K ). 

Hence the argument in the previous paragraph implies that the thermal evolu- 
tion of the primordial gas cloud is determined by the "shock diagram", unless the 
system becomes gravitationally unstable or some processes originating from external 
radiation fields become important. The dynamical evolution of the system is roughly 
determined by the sound crossing time t sc , free fall time iff, and t coo i- In the case 
that t sc is the shortest time scale, the evolution of the system is nearly isobaric, and 
in the case that i coo i is the shortest, the evolution is approximately isochoric. In 
both cases, we can apply the shock diagrams on the thermal evolution of the system. 
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Fig. 5. Same as Fig. 1, except that the cooling time is estimated by assuming the isochoric evolution 



In the third case, when tg is the shortest, the evolution is different from both of the 
above two cases. However, in a realistic situation, such as the shocks in dynamically 
collapsing primordial gas cloud s, at first the isobaric case is realized generally, as 
discussed by Yamada and Nishi. Results of numerical simulations also imply that 
the shock diagrams can be applied in this situation. In fact, in dynamically collaps- 
ing primordial gas clouds, hydrogen molecules are formed up to j/h 2 ~ lCP 3 ,™'^ 
though the gas is not always isobaric because of a non-steady state of the shock 
wave. 

The gravitational instability of sheets may change the thermal evolution of the 
primordial gas from the predicted path in the shock diagram. The shocked and 
compressed layer may fragment into small pieces with a roughly local free-fall time 
which is proportional to p~ 1 / 2 (Ref. |6|)). In steady shock waves, Yamada and NishiEf 
discussed this problem by comparing the time scale of fragmentation and the cooling 
time. 

The external UV radiation field, which exists at the early epoch of the universe, 
may affect the thermal evolution of pregalactic gas. The reaction rates of photoion- 
ization process and photodissociation process have different dependence on density 
from the other reactions induced by collision. As a result, the above arguments based 
on Figs. 1 and 5 are changed significantly. However, the external UV radiation field 
does not exist, as long as we are interested in the first objects in the universe. 

Consequently, once the primordial gas cloud is heated above 10 4 K and there 
is no external radiation field, significant amounts of hydrogen molecules are formed 
(yn 2 ~ 10 -3 ), and the temperature drops much below 10 4 K, until some other time 
scales become comparable to the thermal and chemical time scales. 



of the cloud. 
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Appendix A 

Reactions 

We include the following reactions to obtain the dynamical evolution of the 
postshock layer in §[|. The reaction rates are given in the listed reference in the 
Table |. 



Table I. Reaction rates 



Reactions 



Reference 



H+ + e + 7 
H + e ^BT+ 7 
H~+ H -^H 2 +e 
3H ^H 2 + H 
H 2 + H -^3H 
2H + H 2 ^2H 2 
2H 2 ->2H + H 2 
H + e -^H++2e 
2H — >H + H++ e 
H + H+->H++ 7 
H++ H ^H 2 + H+ 
H 2 + H+-^H++ H 
H-+ H+^H++ e 



H++e 
H++ H"- 



^2H 
H + H 2 
H"+ e -^H + 2e 
H"+ H -^2H + e 
HT+ H+^2H 



Spitzer 1956 1 
de Jong 1972i 
Beiniek 1980E> 
Palla, Salpeter fc^tahler 1983EJ 
SKEf ,_. 
Palla, Salpeter &jJ3tahler 1983EJ 
SKEf 
Lotz 1968 LJ |-, 
Palla, Salpeter & Stahler 1983EJ 
Ramaker & Peek 1976 B 1 
Karpas, Anicich & Huntress ,U|79Eil 



18) 



Prasad & Huntress 19§ n l is l ) 

Poularert et.al. 1978 
Mitchell & Deveau 1983 
Prasad & Huntress, 1980 

Duley 1984B 
Izotov & Kolensnik 1984E 
Duley 19840 ) 



1! ) 



Appendix B 

Steady Shock 

In this appendix, the basic equations of the steady shock are presented. The 
numerical results presented in §|| are based on the following equations. SK solved 
these equations coupled with the chemical reactions. One of the most important 
features of the steady shock is the approximate isobaricity of the postshock layer. 
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B.l. Basic equations 

The jump conditions at the shock front are 

Pivi = P2V2, (B-l) 
Pivf +pi = P2V2 + P2, (B-2) 

(B-3) 

These equations constitute the conservation law of the mass, momentum and energy 
across the shock front. The quantities whose suffixes are 1 and 2 denote the hydro- 
dynamical variable in the preshock region and the postshock region, respectively. In 
case that the shock is strong, i.e. the sound velocity in the preshock region is much 
smaller than the bulk velocity u\, we obtain from Eqs. (BT)~(B-3), 

pz/pi = V1/V3 = 4, (B-4) 
3 

kT = —pm p v 2 , (B-5) 

where p and m p are the mean molecular weight of the gas and proton mass. 

Behind the shock front, the assumption that the flow is steady allows us to write 
the conservation of mass, momentum and the energy as 

pv = p 2 v 2 , (B-6) 
pv 2 +p = p 2 v\ +p 2 , (B-7) 

§ = <r- A) + [<„ + *>]!§, (B.8) 

where r, A and E denote the heating and cooling rates, and the internal energy per 
unit volume, and d/dt is the Lagrange derivative. The hydro dynamical variables with 
no suffix denote quantities at a general point behind the shock front. In addition, 
we require the equation of state p = nkT. Then the equations are closed, and we 
can solve the flow equations. 

B.2. Isobaricity of the postshock layer 

The assumption that the flow is steady indicates the approximate isobaricity of 
the postshock region (Ref. |30|)). From Eqs. (B-7) and (B-8) and the equation of 
state, p = pkT/(m p p), the relative compression factor x = p/p2 is given by 

JlL x 2 _ (1 + _P2 )x + y 2 = Q (B . g) 

rripp p2V$ 
In case the shock is strong, the above equation becomes 

kT 



m p p 

Its solution exists if the condition 

n. n kT 



x 2 - Avix + v$ = 0. (B-10) 
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is satisfied, and this condition holds unless the heating rate is not too high to raise 
the temperature by a factor of 4/3 in the postshock flow. 
The solution is expressed as 



If we define T2 as 

kT2 = 3/j,m p v%, (BT3) 
the pressure p can be written as follows: 

* = V»{ 1 + f-z£)- (B14) 

Thus, p/p2 can be determined by T/T2 only, and the dependence is weak. 

After the gas has cooled down to T <C T2, p has the asymptotic value, 4/3p2> 
which is not much different from p2- As a result, the postshock region is almost 
isobaric for a wide range of temperatures. 
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